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TECHNICAL MEMOEANLUM NO. 359 



NOTE OH THE CALCULATION .OJ 30TJNDAEY LATEES* 

By L. frandtl 



STJMilAET 



The properties of the solutions of the hydro dynanic 
equations of viscous fluid "by " "boundary-layer oaission" 
are discussed. A nethod is indicated for the numerical 
determination of the solution for a known initial profile 
u(x ot y) and pressure distribution p(x) within the region. 

§ 

1. INTEODUCTION 



In flows with snail viscosity (lar/*e Eoynolds nunber) 
the field of viscosity is frequently restricted to a snail 
zone in the neighborhood of the body surface. To facili- 
tate the calculation in this case, the following assunp- 
tions nay be nade: 

1) In tho frictional zone, torced "boundary layer," 

only the strongest terris of the viscosity 
force (jiAu are taken into account (thus, for 
instance, on the x component of this force 
HAu - when the principal direction of flow 
parallel to the wall is chosen as x direc- 
' tlon, and that at right angle to it as y di- 

3 s u 33u 
rection - \x and v are neglected 

ox 3 3z 3 

relative to ' jj, . 

3y B 

2) Outside of the boundary layer the flow is treated 

as frictionless and, because of the snail ox- 
tent of the boundary layer et right angle to 
the wall, the pressure in this layer is equated 
to that which the related fricticnless flow 



*"Zur Eerechnung der G-renzschichten. Zoitschrift fur 
angowandto Mathbnatik und Mochanlk, vol. 18, no. 1, Job. 
1938, pp. ^-83. (Special Eeprint) 



2 



ITACA Technical .Mecorandum ■ No. 359 



uould furnish on the wall.* -Tho trifling 
notion of tho prsBBure with the wall distance 
within the "boundary layer is disregarded, 
hence the pressure considered ' as a given func- 
tion of tho coordinates (x end 2}. , 

Per an average speed u_ of tho ov.ter flow and a 
flow lougth I, the calculation gives for tho "boundary 
layor thickness" 6 tho order of magnitude 

I / — ^— = I / ,/lT" (u = |i/p = kinccatic viscosity). 
,/ u 0 l 

jTor tho two-dimensional caso with x = arc length 
neasurod along tho wall in flow direction ar.d y = verti- 
cal distance on the wall, tho equations road: 

du du . du 1 du d a u N 

— + u r- + v + — ~ = -0 3— r- (1) 
at 3x cy p dz 9y3 



8u 9v . . 

— + — = 0 (2) 
ox oy 

Cf ger.oral interest is the stationary caso ( ^ = 0 N ) . 

^3t / 

about v/hich nuch as boon writton (see reference l) , al- 
though thcro still is no satisfactory solution for one 
important aspoct; lamely, tho further development of a 
given velocity jrrofilo by given pressure distribution. 
This problem forms the snbjoct cf the prosont discussion. 
In mathematical language it is cxproseod as follows: 

Given: the velocity profile u(x Q ,y? for tho Initial 
section x = x 0 , and, in addition, 

^ f£ = T,f(x) for - x ^ x 
p dx • 0 

To find:. u(x,y) for x > x n . 



♦More ;prociselj , the f rict ionless flow should bo takon on a 
body thickened up by the displacement thickness 

6* = — / (u x - u) dy. : 
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2. PEOPEBTIES OS THE SCATIOEABT BOUITOABY LATXB SLOW 



The known solutions contribute nothing to this problem 
because they contain noroly expansions in povors of x 
fron an initial section, vaoroin u = 0 or u = const 
(stagnation-point flow and.' flow past an edgo with aero 
edge anglo). On the other hand reforonce nay bo node to 
tho fact that - discounting the caao of brook-away of 
flow - tho oiistonco of tho solution is -physically evident : 
a boundary layor profilo croated during any procoding 
poriod nust, so long as no separation occurs, continuo to 
dovclop in Bono way or other. Tho mathematical answor to 
this question, which I already hold at tho first mathe- 
matical congress at Heidolborg in 1904, is as follows:* 

By t.io uso of oquation (2), u^i- + v ~ in aquation (l) 

nay bo rowrltton as: v -i^i — u ^21 = - u 3 -S- f ~\, honco, 

oy dy oy v. u /* 



5- " ^ ^ (f(x) -0) *7 + 9(x) (3) 

* is .tho slopo of the stroanlino toward tho x axis. 

If tho axis itsolf is tho wall, thon ~ = 0 for y = y 0 = 0, 

i.e., o ( x) * 0. As a secondary product, it is soon that 
in oquation (-3) the caso of a (slightly) wavy wall is also 
includod, which is different fron tho x axis. It is, 
then, 

• : T&.-(i) ■»(-) 

* « O 



•In tho toxt of tho sposch at tho III Int. Math. Cong., 
Leipzig, 1905, and in tho new odition by PranaVtl-Bot z on- 
titlod: T"ior Abhandlungon zur Eydro&ynanik und Aerodynanik, 
GBttingon, 1937,.. p. 487, it roads: if, as usual, dp/dx 
is given, and, in addition, the curve u for tho initial 
sootion, cvory such problon can bo numerically dealt with, 
by obtaining through quadratures fron oach u the rolativo 
da/ ox; in. this nau'ner it "is always' possible to novo a stop 
forward in tho x diroction with tho aid of one of tho 
known approxinato nothods. A difficulty exists of course 
in different singularities occurring on tho fixed boundary. 
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which leads to 

7 0 U) = I* <p(i) dx 

.is oquation for the wall. This secondary result can also 
bo expressed so that, if u(x,y) is a solution of tho 
syston (l) , (2), u(x,y - y Q (x)) itsolf is a solution. 
Haturally, this can also "bo provod direct fron (l) , (2). 

The significance of aquation (3) for our problon, 
howoTor, rocts on the fact that 6v/fl;' and hence 3u/3x 
can "oo obtained by ncar.s of equation (2). Per <p(x) = 0, 
it is 

• *f5>*j (4) 

Tho right-hand .sido of this equation r.erely contains dif- 
ferentiation and intogration with respect to y, and is 
thorbfero ar:onablc to couploto calculation for ar. initial 
profile u = u(x 0 ,y), fcr which it yields 3u/3x and with 
it tho possibility of continuing with ropsoct to x ( pro- 
vided that 4!r is restricted. It gives a new u for 
ox 

x 0 + ax: 

u(x 0 + ax, y) 

Tho r.othod can bo repeated for this new profile and yields 
a u'x + 2 a x, y) otc. f thus proving tho oxistor.ee of 
tho solution for all zones wherein u renains positive* 
On the other hand, it is readily apparent that tho nothod 
fails if u passes anywhbre through zero, whero tho in- 
tegrand becomes quadratic ally infinite. 17ow, however, 
u = 0 along tho wall, that is, for cp(x) ■= 0 along axis, 
x, because of tho adhorenco of tho fluid to the wall. Tho 
condition is ameliorated becauoo v tocor.os at tho sane 
tine quadrat i sally zero.* I7ow, a glance at equation (3) 



= u(r. 0 ,y) + |^(x 0 ,y) Ax 



•With u, |^ =0 also, for y = 0; but v = 
bocauso of equation (2)-. 



.0 
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indicates that tho function ''is regular only when tho 

lift o grand '1b everywhere finite This - stipulates that, 

a . 

on the wall, tho expression f(x> -.-^-j quadratic ally 

a ^ 3 

diaappoars, i.o., that - t ^ and 



X=X, 



eiT ° 9 - K^/C^T) T 

X-X ( 



Tho evaluation of tho "lin" then 

for which ocuation (4) ■ gives tho Uniting valuo of 

.-■ -I o 

this r.cr.-ia that 

1) The second and tho third dif f orcntial quotients of 

tho initial profile for 7=0 cannot he arbi- 
trary, hut rather nust "be = f( x 0 ) m& 0. 

2) The fourth differential quotient at y «= 0 is do- 

cisivo for th 3 further dovolopncnt of tho pro- 
file. 

Por a contor.platcd second stop, thoroforo, tho initial 
profile u(x D ,y) +^^^1 a:: alroady contains — She 

i=x 0 3 y 

now "lin" consideration shows that tho fifth and the sixth 
differential quotients aro again hound and tho success of 
the socond sto; in wall proxinity doponds on the Boventh 
differential quotiont of tho initial profile at y = 0J 
How it nust he considered how, for instance, tho seventh 
differential quotient is to "bo deterrinod on a tabulated 
function, and at tho one end cf tho table of functions at 
that. And all that is accomplished hore are but ttco .short 
steps. 

On top of that tho bonds.' It is soon fron oquation 
(3) or (4) that a violation of the- bonds loads to a singu- 
lar bohavior of tho solution at points x =' x 0 and y = 0. 

Two questions arise: ■ • 



l) How is tho vicious action of the oquation syston 
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(1), (2) in the vicinity of places with u = 0 
to do understood 

2) Of what typo are the s inguinrit ios which occur on 
a regular, but tho bonds violating, initial 
profile ? 



3. P". OCZD'JUE Or SOLUTION AC PLACES VflTE u < 0 



At this instance it is first to ba stated. that tho 
ITavior-Stokes equations r.ro closoiy rclatod to tho bihar- 
nonic equations (LAT = 0) , whenco tho solutions within the 
zone aro always analytic. Tho introduction of the bound- 
ary layer onissione, hcv.-ovor , represents a character- 
changing intervention. Tho structure of equation (3) it- 
self indicatos the prosonco of characteristics (nanely, 
the straight linos parallel tc axis y) along which dis- 
continuities occur, if f(x) or 3>(x) havo discontinui- 
ties. This points therofcro to a paracolic character. 
3y a contact transf crnation by noans of which the stroaa 
function \Jr is introiucod as an independent variable,* 
tho case bococes clearer. Posting 

u = ^ and v = oauation (2) is identically satis- 

fy ox 

fiod. liriting the now coordinates | , ti in place cf x,y 
in oauation ( l) , and posting £ = x , T] = y gives 



•In the spring of 1914 I was occupic:! with the flow through 
narrow channels and. found, when applying the usual boundary 
layor equations, soao difficulty in fornuliting tho Unit- 
ing , conditions u = v = 0 on tho socond e*'.go, I thorcforo 
attenptod to introduce ths strsan function as an independent 
variable, since the two ed^ta now siiaply read V = 0 and 
\|/ = Q = court. Equation (5) was subsequently rediscovered 
by Von r-iiscs and reported at the Eissingor neoting of the 
Society for Applied Mathematics and Ilechanics in the fall 
of 1927 (Z.i .a.H.K. , vcl. 7, 1327, p. 429). Since I failod 
at that tino to publish ny own calculation rothod. Ton Icises 
has thcroforo the usual priority. Jor tho rest, tho rcador 
is referred to ny rcnarks incidental to Von iliscs ' report 
in this periodical, vol. 8, 1923, pp. 249 and, particulary 
section 2, 250. .: 
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"total pressure" 
"being disregarded) introduced herein and x substi- 



f or 



| t and posting u p = |i gives 

a§ -.. .. a^g 



I 

(•*• 



o x 



= H u 



3ya 



then 



(5) 



whore u = ,/3(g - p)(z))/p. Squation (5) is then a differ- 
ential equation which is closely rolatod to that of hoat 
conduct! or. in place of bar length and x instoad of 
tino) , nr>t\ of tho parabolic typo of which is Icnown that on 
it a propitious direction with analytical disappearance of 
tonporaturo dif f oro:-.cos , etc. ( incr oas in^ tine periods) 
exist, but at the sc.no tine also a deleterious direction 
with occurrence of singularities (extrapolation to states 
which with respect to tine are before tho givon initial 
distribution).* On oquation (5), sinilar to tho problon of 
hoat conduction, the dovolopnont with respect to increas- 
ing values of x is propitious so long as u- is positive 
But whore ]i is negative, as bohind tho broak-away point, 
tho dolotoriouB stato begins. Ph;'sically this noans that 
tho wall-adjacont fluid particlos bohinr" tho broak-away 
point move fron groator toward snnllor z, whonco thqir 
volocity distribution doponds upon tho phenonona'n taking 
placo at theco 'greater x values. Attempts to defino 



*"Eow should tho distribution havo looked before a ninuto 
(or hour) in ordor that tho nonontary nonuniform tonpora- 
turo -distribution originates from it" (ovidontly froquontly 
an unsolvable problon J) . 
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th.cn fron tho states, at snallor x Is, at host on a par 
with guossinj ganos. In other words: the corroct subse- 
quent dovclopnont of the flow conforning to obsorvationa 
aft of tho "broak-away point™ ia prcdlcatod on tho caroful 
soloctlon of cur function f(x) to fit tho desired procosa. 

As concerns the behavior of tho boundary layor on tho 
wall It30lf in tho "nomal rogion" (u > C for y > 0) , 
the analytical continuation of u In tho rjgion of no.ja- 
tivo y loads to negative values of u, if, as stipulated, 

~ ia -positive for y = 0. Tho tondoncy of tho solution 
oy 

toward singularities along the wall is thus explained by 
tho fact that tho analytic--! continuation through tho wall 
loads into a "dotrinontal zonoM To tho oxtont that power 
scries with respect to y are available for separate solu- 
tions, tho extrapolation car. be achieved. ■ It showB in tho 
zou of negative y a continuous ascent of tho function 
values, such as y o a I"* or slpilar functions. 

a. solu!pio:: 3Y yioiatios op the n 3o:rDS B 



Tor a Ion.; tine I unsuccossf v.lly sug-ostod the problon 
citod at tho ond of section 2 as a thesis for a degree of 
doctor until ia 1923 Dr. S. Goldstein, Ganbridgo (England) 
as scientific guoBt of tl;o Institute undertook this problon 
ani explored it thoroughly. Tho singularities inducod by 
an initial profile violating "bonds" are of algebraic na- 
ture It is found that every . initial profile dovolopablo 
in powor scrios for y furnishes for a given dp/dx as 
whole "function a solution which appears as powor scrios for 

^x - x 0 ; tho powor coefficients thcrsselvos boing functions 

y 

of 1 - ■ defined by differential o^uations cf tho 
-A - x Q 

third order. (Details nay be found in Goldstein's roport 
ontitlod "Concerning sono solutions of tho boundary layor 
oquations In hydro lynanico , " Proc. Canbr. Philos. Sec, 
vol. 25, I, 1330, p. 1.) Tor tho problon In hand this 
nothod of solution does not appear tc be Innediatoly appli- 
cable. And cvon if It cculd bo continued, the gain would 
not bo nuch because 3n functions would havo to bo fornulatod 
in ordor to reach the nth power of x - x 0 and that In- 
volves a furthor lncroaso in tho anount cf papor work with 
ovory power of fa - x Q . Goldstoin. ^figured with throo torns 
and ovon then required vory conplicatod dovolopnonts in or- 
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der to establish the connection with the given function 
Tallies.* He thus reached only the linear tern, whence the 
"stop -Ax' ■ would have' to remain . very, anall so as not to 
introduce inadmissible orrora. The new profile would then 
have to to represented again as power sories of y in 
order to prepare tho next stop. The whole procoss would 
thoreforo bocono extremely laborious and hardly profitable, 

Thoro is thoreforo thq greatest intorost for continu- 
ation with integral powors of x. And tho subsequent argu- 
ment nakos tho prospects for -a practical way toward this. 
ond t appear hopoful. If it involved tho continuation of 
a solution already achiovod by a provious method f or a 
completely rogular c our so of tho prossuro gradient, tho 
collective bonis should, by rigorously effoctod calcula- 
tion for tho prosoding eoction air only bo satisfied in tho 
o:.d profilo which is to servo as initial profile of tho 
now section. They aro net exactly such in reality, for 
the scrios expansion in x is either stopped aftor a com- 
paratively lev/ torn or ols*"o numerical inaccuracies have 
occurrod. So, ordinarily, it will corcly becono a mattor 
of roostabliehing tho bonding conditions by vory slight 
changes in tho givon initial profile, after which tho con- 
tinuation in integral powors cf x nay proceed. 

Tho question, hov eventual errors undo by this shift- 
iiir- of tho profile affoct tho subsequent stops, is in order, 
this will involve mostly such departures from the unknown 
oxact for::, now winding positive, thon nogr.tivo, around 
the enact curvo. According to tho theory of heat conduc- 
tion, such short-wove fluctuations in tho tenporaturo dis- 
tribution dio out relatively quietly. Honco only tho prox- 
imity cf tho v?„lue u = 0 noods spocinl study, and for 
this very reason I suggostod the case of a snail disturbance 
€ of oscillating courso superimposed on a velocity profilo 
of the simple form u = a y. Stipulating, for reasons of 
sinplif icntion, that e shall bo of tho forn c = f(x) cp(y), 
it follows fron systems (l),(2), aftor linearization, that 
it is necessary to put f(x) = o _ax . Tor <p(y) tho differ- 
ential oquation thon gives 

a a y qp*(y) + vo'"' (y) = 0 
which is colvod ty 
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(This solution was supplied by my co-worker, Dr. H. GSrtler. 
He also adduced the proof that a solution for f = e +ax 
satisfying one of the limiting conditions does not exist.) 
Of these expressions, however, only that with A satisfies 
the Uniting conditions cp(y) = 0 and <p M (y) = 0; hence 
the physical solution in question is restricted to this 
term. The power expansion of y gives a term with y, 
another with y 4 , one with y 7 , etc., as corresponds to 
the "bonds for the case f(x) = 0. As for the rest, the 
function oscillates with decreasing amplitude and decreas- 
ing wave length for increasing y.* This oohavior of the 
solution may therefore he cited as proof that, "by any small 

errors of (g^^ la tho Tel °city profile omployod for 

continuing the solution, no aftereffects oxecoding tho 

amount of tlioso 

is positivo. 



errors are oxpoctod so long as 



V oyy y =o 



5. TEE TOW M3TZ0D 

In ordor to gain r. olonr plan for the final mothod, 
the system of tho "bonds is first oxplorod "by moans of an 
expansion of u in powers of y "by insertion in tho sys- 
tem (l), (s). Suppose that 

a s y 3 av.y k 

u = -iy + —^r +• •+ ~tT~ + ••• 



whero a 1 ...aj E roproaont functions of x. Differentia- 
tions with rospoct to x "being donotod "by dashos (qj', 

etc.); — = y f(x). It affords in ordor tho following 
. p dx 

interrelations, if oach a n with n > 1, whorever it oc- 
curs on tho right-hand side, is replaced "by its vcluos 
given sevoral lines abovo it: 



•Tho asymptotic fcrnula for J ±j./3 6* T °s 1 disappoaranco 

constant a = ^ _ nnn r , wfV vo longth 7\ proportional 

ay(\/2n) 

. 1 
to — *> • 
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a x (free), a 3 = f (givon) , a 3 = 0 (given), 

a* - ^ oiGi 1 (f reo) , a B ■ — ajf 1 (given x a ) , a a = — ff 1 (given) , 
a ? B ^i^ 4 &:ia ai ' * r , 'a 1 'ai u J -(froo) 

If tho ecrios is stopped at this point, the throe freo 
values a x , a 4 , a 7 can be obtained from throe linear equa- 
tions, by proceeding froa tho u values rolatod to throe 
fixed valuos of y. Sho first to tho sovonth powers of tho 
three y values are computed boforohand and usod again for 
every stop. Posting a x in a 4 givos aj 1 , a x and a x * 
in a 7 c 1 * tfith those values 

a 4 1 = i (a 1 a 1 M + a x ' a ) can still bo conputod. 

I2:o v/holo method can bo refined by subsequent conpu- 
tation of tho torr.s of eighth and ninth pov/ers -.rith 

a 8 =_L|lO a x a f" - 13 a 1 a 1 » f + 9(a 1 a{»+ a^a) f 
and 

a b = ~ [40 a.ff 1 ' - 16 a, f ,s 1 

for the tlirco y valuos, and thon by calculating ag-.in 
a lt a 4 , and a 7 froc tho givon throo u valuos (to bo 
ropoated, if necessary, by iteration). Tho cost appropriate 
choico of the throe y values is subjoct to a spocial 
study. It will likely bo found appropriate to procood by 
first obtaining a first approximation of a x froc. an ordi- 
nate 7jj then of a 4 by r.oans of j r x and y , and an 
improvement of o lt lastly a T with three y values, an 
inprovoaont of a 4 and a socond improvement of a t . By 
this method tho second. Btop can bo extondod to a 6 , tho 
third to a 9 , tho proco-Iing approximation of aj and a 4 , 
respectively, being each tine posted in- tho higher torms, 

ITond icons icnal coordinates nay bo chosen for the 
numerical calculations, by which, in place of v tho unity 
(X = x/l, 

T = 7 </u 0 /u I ', U = u/u b ) 

substitutes. 
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! This nothod does not furnish a 7 '; hence it cannot ho 
continued without contact with another nothod which sup- 
plio's tho dovolopnont of tfho voloclty profilo on tho insido 
of tho fluid; For it, howovor, tho rulo given in oquation 
(4) or olao that of oq.ur.tion (5), which probably is ovon 
noro suitahlo, is availablo. This noroly noods to bo con- 
plotod by a transition between y and . 

Eo cause of 

u =Ji ( g - P ) = || 

wo havo 



Posting heroin 



to bo inprovod, if nocossary, analogous to tho Rungo-Kutta 
nothod) , tho y value related to cvory \j/ or u is known 
up to tho unknown valuo 0 (s; 0 + Ax). This donotos a par-- 
allol shift of tho velocity profilo obtainod and can bo 
so dotcrr.incd that tho profilo continuously Joins tho wall- 
adjacont profilo obtained by tho other rethod. Provision 
for a wide ovorlap of both profiles nust bo nado, which 
thus affords tho r.oans for tho dotcrninat ion of tho coef- 
ficient a 7 for tho next step. Conversely, tho connec- 
tion for the continuation on tho inside is also important, 

since tho nunorical determination cf g on tho odgo of 

tho region would each tine bo accompanied by a loss in 
width unless tho course of g in the edge strip is in- 
cluded. 

ITaturally, tho stop of x 0 + ax can be oxocutod also 
aftorward in the wall-ad jaccnt. strip backward toward x 0 
according to the old nothod and a correction of the socond 
ordor dorivod fron the deviations' (if necessary also four 
stops bach and forth, so aa to lowor tho orror of tho stop 
to tho fifth ordor, as in P.ungo-Kut ta ' s nothod, which thon 
naturally onablos the choico of groatcr steps Ax). 

Tho nothod, ovon as doscribod, is not oxactly sinplo, 
but it should givo noro reliable results. Its practical 



y = 5(x) + /""' d v/ u 
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approval Is hol:l in aoofanco ponding tho findings of tho 
. £uisor.Vilhgln Institute for Plow Sosoarch. It is expect- 
ed that tho continuation can bo carried on to tho "break- 
away point and probably a littlo boyond. A study of tho 
conditions on tho break-away point (a x = 0) shows that 

a x ' assumes tho forn sinco a 4 along with a 1 bo- 

conos zero. But continuation is possiblo if a x 1 and a 4 1 
which hero boconos = a 1 ,a /i>t' is oxtrapolatod as far as 
that place. Iron tho break-away point x = x^, it prob- 
ably will bo posciblo to penotrato a littlo farther with 
an oxpansion in powers of x - x^ and y t tho coeffi- 
cients of which aro given by tho continuous connection 
with tho region x< x_^, whore, by that tir.o the courso 
of f(x) will havo to be cheson cn the pronisa that a 
dond-air rogion shall occur. 



Translation by J. Vanicr, 
national Advisory Connittoc 
for Aoronautics. 
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